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Abstract 

We consider the T-Stein equation X = AX T B + C, where the operator (-) T 
denotes the transpose (T) of a matrix. In the first part of this paper, we 
analyze necessary and sufficient conditions for the existence and uniqueness of 
the solution X. In the second part, a numerical algorithm for solving T-Stein 
equation is given under the solvability conditions. 
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1. Introduction 

The purpose of this article is to analyze the solvability conditions of the 
T-Stein equation 

X = AX T B + C, (1) 
where A, B, C E i£_mxn are known matrices, X G c"ix« j g a ma ^ r j x i Q h, c 
determined. Our interest in the T-Stein equation originates from the study 
of completely integrable mechanical systems, that is, the analysis of the T- 
Sylvester equation 

AX + X T B = C, (2) 

where A, B, X are matrices in R mxn [l|,|2|. By means of the generalized inverses 
or QZ decomposition Q , the solvability conditions of @ are studies in 0, 0, H| ■ 
Suppose that (A, B T ) is regular, that is, aA + bB T is invertible for some scalars 
a and b. The T-Sylvester equation ([2]) can be written as 

{aA + bB T )X + X T (aB + bA T ) = aC + bC T . (3) 
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Premultiplying both sides of © by (a A + bB T ) 1 , we have 

X + UX T V = W, (4) 

where U = (aA + bB T )-\ V = aB + bA T and W = (aA + bB T )- 1 (aC + bC T ). 
This is of the form (fTJ). In other words, numerical approaches for solving 
can be obtained by transforming © into the form of ([T]), and then applying 
some iterative methods to ([T]) for the solution 0, H, i| . With this in mind, in 
this note we are interested in the study of T-Stein matrix equation (TT]). 

Observe that if X is a solution of ((T|), then X also satisfies the following 
Stein equation 

X = AB T XA T B + C + AC T B. (5) 

It follows that once there exists an unique solution to JS]), we can usually solve ([T]) 
via the conventional Smith-type iterative methods given in 7, 8] to ©. How- 
ever, even if ([T]) is uniquely solvable, it does not imply ([5]) is uniquely solvable. 
For example, if b = 1 and a = — 1, the scalar equation x = —x + c has an 
unique solution x = — , but no more information can be obtained from the Stein 

equation x — x + c — c. In our work, we formulate the necessary and suffi- 
cient conditions for the existence of the solution of ([T]) directly by means of the 
spectrum analysis in Section 2. We study the capability of the Smith- type iter- 
ative algorithms for solving ([1} in Section 3 and concluding remarks are given 
in Section 4. 



2. Solvability conditions of the Matrix Equation X — AX T B + C 

In order to formalize our discussion, let the notations A®B be the Kronecker 
product of matrices A and B, I n be the n x n identity matrix, and <j{A) be the 
set of all eigenvalues of A. With the Kronecker product, (P) can be written as 
the enlarged linear system 

(J ron - (B T ® A)V)vcc(X) = vec(C), (6) 

where vec(X) stacks the columns of X into a column vector and V is the Kro- 
necker permutation matrix [9( which maps vec(X) into vec(X T ), i.e., V — 
Y^, e j e J ® e^ej, where ej denotes the i-th column of the mn x mn iden- 

l<~i,j<mn 

tity matrix I mn . Due to the specific structure of V, it has been shown in [ToL 
Corollary 4.3.10] that 

V T {B T (g> A)V = A® B T . 

It then follows that 

{{B T ® A)Vf = (B T ® A)VV T {A ® B T ) = B T A ® AB T , (7) 

since V 2 = I mn and 7? = T T . Note that cr(yl T B) = a(B T A) = a(AB T ). 
By ([7]) and the property of the Kronecker product [llj, Theorem 4.8], we know 
that cr(((£ T ® ^)7 5 ) 2 ) = {AjAj|Ai, Aj € a(A T B) = {A x , . . . , A„} , 1 < i,j < n}. 
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That is, the eigenvalues of (B T ® A)V is related to the square roots of the 
eigenvalues of a(A T B), but from (|7|), no more information can be used to decide 
the positivity or non-negativity of the eigenvalues of (B T ® A)V . A question 
immediately arises as to whether it is possible to obtain the explicit expression 
of the eigenvalues of (B T ® A)V, provided the eigenvalues of A T B are given. 
In the following two Lemmas, we first show that the generalized inverse of an 
upper triangular matrix is still upper triangular and then apply it to discuss the 
eigenvalues of (B T ® A)V. 

Lemma 2.1. Suppose that A £E C mx ™ j s an upper triangular matrix having full 
rank. Then the generalized inverse of A is an upper triangular matrix. 



Proof. Let A 



A 1 



, if to > n, or A 



A 



, if m < n. Here, A 



is an upper triangular matrix with rank(A) = rank(A). Then, the generalized 



inverse A> of A is an 



matrix satisfies A' 



A- 







, if m > n, or 



A< = A 



l-T 



ifm < n. This concludes that A^ is upper triangular. 



1. There exist unitary matrices P G C" 1 *™ an d 



such that Ua 



PAQ and Ub '■= Q H B T P H are two upper triangular matrices. 

2. (B T ® A)V = {Q® P H )(U A ® U B )V(Q H ® P) 

3. (j{{B T ®A)V) = {Ai, : v Xa7A,.A, e a(A T B) = {Ar, . . . , A n } , 1 < i < j < n}. 
Here, *Jz denotes the principal square root of a complex number z. 

Proof. Assume first that B is full rank and PAB T P H = T\ is the upper schur 
form of AB T . This defines the unitary matrix P and one can consider the QR 
and QL decomposition of matrices B T P H and (PA) H , respectively, such that 

QH B TpH = T ^ f orm > n; 
Q H {PA) H = f 2 T , form < 7i, 

where Q and Q are unitary and Ti and Ti are upper triangular. By Lemma l2.ll 
there exist two generalized inverses T 2 and T\ of matrices T2 and T2, respec- 
tively, such that T%r\ = I n and T^Ti = I m . In turn, one shall have the following 
upper triangular matrices: 

U A :=PAQ = TiT%, U B ~ Q H B T P H = T 2 , for m > n, 
U A :=PAQ = T 2 , U B := Q H B T P H = f^T 1 , for to < n. 

If B is not full rank, one then considers the SVD decomposition of B: 

B = Udia,g(a ir -- .oy, (),••• ,0)^ T , 
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where U G c mxs arLC ] y £ i^sxn arc un jtary, s = min(m, n), and r = rank(B). 
This then defines a sequence of nonsingular matrices 



B *--= B + \ E ^ T < 



i=r+l 



for any positive integer k. It follows from the discussion above that for any 

xm and Q k G C" xn 
fc B^,P" are two upper triangular 



P 



< 



positive integer fc, there exist unitary matrices Pk G 
such that U Ak := P fe AQ fc and P Bfc := <^h u * v h 
matrices. Since ||P fc || = ||Q fc || = 1, ||?7aJ|f = \\M\f, and \\U B 
\\B\\p + ~tt~ ; h follows from the Bolzano -Weierstrass Property [12J that there 
exists a convergent subsequence of (Qk, Pk, Pfe) such that its limit is equal to 
(Q, P, P) for some unitary matrices Q and P, which completes the proof of 
part 1. 

From part 1, there exist two unitary matrices P and Q giving rise to 



(B 



)A)V = (Q<E>P H )(U B <E>U A )(P<E)Q H )V 
= (Q®P H )(U A ®U B )V(Q H ®P). 



Thus part 2 holds. 

Let the diagonal entries of Ua and Ub be denoted by {an} and {bjj}, re- 
spectively. Then, (Ua ®Ub) is an upper triangular matrix with given diagonal 
entries, specified by an and bjj. Here, one can assume without lose of generality 
that Ua and Ub are n x n matrices, i.e., m = n. After multiplying (Ua ®Ub) 
with V from the right, the position of the entry aubjj is changed to be in the 
j + n(i — l)-th row and the i + n(j — l)-th column of the matrix (Ua (8 Ub)V. 
They are then reshuffled by a sequence of permutation matrices to form a block 
upper triangular matrix with diagonal entries arranged in the following order 



aub 



By 



an 6 22 
a 2 2&n 

it can be seen that 



, a 2 2h 



Oll^rm 

a nn b n 
q n — 1 _ n —\b nn 

^nn^n—l.n—l 



a 2 2^>33 



(8) 



cr((P T <g> A)P) = {«„/)„. : v u, , a ,,!>,, !>,,. 1 < i,j < n| 
= jA^i^AiAj, 1 < i,j < ?ij 

where Xi = aubu G <r(A T B) for 1 < i < n. Note that if to ^ rt, the diagonal 
entries given in ([8]) will remain unchanged, except that some zero diagonal 
entries will be added or removed from the diagonal. 



Note that when m = n, the decomposition given in part 1 of Lemma 12.21 is 
called the periodic Schur decomposition The following result, providing the 
unique solvability conditions of ([I]), is an immediate consequence of Lemma l2.2l 
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Theorem 2.1. The T -Stein matrix equation ([T]) is uniquely solvable if and only 
if the following conditions are satisfied: for A G a(A T B) and A 7^ —1, — G" 
o~(A T B); — 1 can be an eigenvalue of the matrix A T B, but must be simple. 

Proof. From (|()]), we know that the T-Stein matrix equation ([T]) is uniquely 
solvable if and only if 

1 o-((B T (g) A)V). (9) 



By Lemma |2~51 if A G <r(A T B), then - £ er(A T B). Otherwise, 1 = J\ - - £ 

A V A 

((_B T (g) On the other hand, if -1 G <r(A T 5) and -1 is not a simple 

eigenvalue, then 1 G o"((£> T <g) ^4)7'). This verifies © and the proof of the 

theorem is complete. 

3. Smith-type iterative methods 

Originally, Smith-type iterative methods [8| are developed to solve the Stein 
equation 

X = AXB + C, A G K mxm , B G R nxn , C G R mxn . 

Our main thrust in this section is to extend the Smith-type iterative methods 
to solve by means of the formula (|S|) and also preserve the convergence 
properties embedded in the original methods. Here we summarize the iterative 
methods as follows. 

• The r-Smith iteration methods: 

r-l 

Xk+i = A l k X k B l k , 

i=0 

where 

X = AC T B + C, A = AB T , A k+1 = A r k , k> 0, 
B Q = A T B, B k+1 =B r k , k > 0. 

It can be seen that if p(A T B) = p(AB T ) < 1, the r-Smith iteration converges 
to the unique solution of (UJ 



X = Y^{AB T y{AC T B + C){A T B)\ 



When r = 2, the r-Smith iteration is also called the Smith accelerative it- 
eration 7] and it has been shown in Q that the Smith accelerative iteration 
performs more effective than any other r-Smith iterations. One possible draw- 
back of the Smith-type iterative methods is that it cannot handle the case when 
the eigenvalue —1 of A T B is simple. Based on the solvable conditions given in 
this work, it is possible to develop a specific technique working on the particular 
case and it is a subject currently under investigation. 
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4. Conclusion 

In this paper, we obtain necessary and sufficient conditions for the existence 
and uniqueness of the solution of the T-Stein equation and find that the solv- 
ability conditions between (QJ and ([5]) are different. When the equation ([T|) is 
consistent, we extend the capacity of Smith-type iteration methods from the 
Stein equation to the T-Stein equation. 
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